Using In-Shoe Inertial Measurement Unit Sensors to Understand Daily-Life Gait Characteristics in Patients With Distal Radius Fractures During 6 Months of Recovery: Cross-Sectional Study.

BACKGROUND: A distal radius fracture (DRF) is a common initial fragility fracture among women in their early postmenopausal period, which is associated with an increased risk of subsequent fractures. Gait assessments are valuable for evaluating fracture risk; inertial measurement units (IMUs) have been widely used to assess gait under free-living conditions. However, little is known about long-term changes in patients with DRF, especially concerning daily-life gait. We hypothesized that, in the long term, the daily-life gait parameters in patients with DRF could enable us to reveal future risk factors for falls and fractures. OBJECTIVE: This study assessed the spatiotemporal characteristics of patients with DRF at 4 weeks and 6 months of recovery. METHODS: We recruited 16 women in their postmenopausal period with DRF as their first fragility fracture (mean age 62.3, SD 7.0 years) and 28 matched healthy controls (mean age 65.6, SD 8.0 years). Daily-life gait assessments and physical assessments, such as hand grip strength (HGS), were performed using an in-shoe IMU sensor. Participants' results were compared with those of the control group, and their recovery was assessed for 6 months after the fracture. RESULTS: In the fracture group, at 4 weeks after DRF, lower foot height in the swing phase (P=.049) and higher variability of stride length (P=.03) were observed, which improved gradually. However, the dorsiflexion angle in the fracture group tended to be lower consistently during 6 months (at 4 weeks: P=.06; during 6 months: P=.07). As for the physical assessments, the fracture group showed lower HGS at all time points (at 4 weeks: P<.001; during 6 months: P=.04), despite significant improvement at 6 months (P<.001). CONCLUSIONS: With an in-shoe IMU sensor, we discovered the recovery of spatiotemporal gait characteristics 6 months after DRF surgery without the participants' awareness. The consistently unchanged dorsiflexion angle in the swing phase and lower HGS could be associated with fracture risk, implying the high clinical importance of appropriate interventions for patients with DRF to prevent future fractures. These results could be applied to a screening tool for evaluating the risk of falls and fractures, which may contribute to constructing a new health care system using wearable devices in the near future.

Arch height flexibility is associated with plantar fascia tension during running.

BACKGROUND: Plantar fascia tension is considered to cause plantar fasciitis, and medial longitudinal arch decrease is believed to be a risk factor for plantar fasciitis. Arch height index (AHI) and arch height flexibility (AHF) are useful indicators for evaluating medial longitudinal arch. However, the relationship between plantar fascia tension during running and these indicators remain unclear. RESEARCH QUESTION: Are the foot characteristics in AHI and AHF that represent medial longitudinal arch related to plantar fascia tension during running? METHODS: Twenty-two male participants enrolled in this study. Foot characteristics required for calculating AHI and AHF were measured using the AHI measurement system. AHI was defined as the height from the floor to the dorsum of the foot divided by the truncated foot length with 10% or 50% load. AHF was defined as the change in arch height from the 10% and 50% loads. Marker trajectories of the foot and force plate data during running were measured using a three-dimensional motion analysis system and a force plate. Based on the measured data, the peak values of the plantar fascia tension were analyzed. Pearson's correlation was used to determine the relationship between foot characteristics and plantar fascia tension. RESULTS: No significant correlation was found between AHI in the 10% load condition and plantar fascia tension (r = -0.36, p = 0.09) or between AHI in the 50% load condition and plantar fascia tension (r = -0.148, p = 0.515). In contrast, a significant moderate positive correlation was observed between AHF and plantar fascia tension (r = 0.568, p < 0.01). SIGNIFICANCE: AHF is a change in arch height between sitting and standing positions, can be easily used to evaluate plantar fascia tension in clinical settings. This study implies that evaluating AHF is a useful tool in considering plantar fascia tension during running.

Gait characteristics in patients with distal radius fracture using an in-shoe inertial measurement system at various gait speeds.

BACKGROUND: Distal radius fractures (DRF) commonly occur in early postmenopausal females as the first fragility fracture. Although the incidence of DRF in this set of patients may be related to a lower ability to control their balance and gait, the detailed gait characteristics of DRF patients have not been examined. RESEARCH QUESTION: Is it possible to identify the physical and gait features of DRF patients using in-shoe inertial measurement unit (IMU) sensors at various gait speeds and to develop a machine learning (ML) algorithm to estimate patients with DRF using gait? METHODS: In this cross-sectional case control study, we recruited 28 postmenopausal females with DRF as their first fragility fracture and 32 age-matched females without a history of fragility fractures. The participants underwent several physical and gait tests. In the gait performance test, the participants walked 16 m with the in-shoe IMU sensor at slower, preferred, and faster speeds. The gait parameters were calculated by the IMU, and we applied the ML technique using the extreme gradient boosting (XGBoost) algorithm to predict the presence of DRF. RESULTS: The fracture group showed lower hand grip strength and lower ability to change gait speed. The difference in gait parameters was mainly observed at faster speeds. The amplitude of the change in the parameters was small in the fracture group. The XGBoost model demonstrated reasonable accuracy in predicting DRFs (area under the curve: 0.740), and the most relevant variable was the stance time at a faster speed. SIGNIFICANCE: Gait analysis using in-shoe IMU sensors at different speeds is useful for evaluating the characteristics of DRFs. The obtained gait parameters allow the prediction of fractures using the XGBoost algorithm.

Foot characteristics of the daily-life gait in postmenopausal females with distal radius fractures: a cross-sectional study.

BACKGROUND: Gait decline in older adults is related to falling risk, some of which contribute to injurious falls requiring medical attention or restriction of activity of daily living. Among injurious falls, distal radius fracture (DRF) is a common initial fragility fracture associated with the subsequent fracture risk in postmenopausal females. The recent invention of an inertial measurement unit (IMU) facilitates the assessment of free-living gait; however, little is known about the daily gait characteristics related to the risk of subsequent fractures. We hypothesized that females with DRF might have early changes in foot kinematics in daily gait. The aim of this study was to evaluate the daily-life gait characteristics related to the risk of falls and fracture. METHODS: In this cross-sectional study, we recruited 27 postmenopausal females with DRF as their first fragility fracture and 28 age-matched females without a history of fragility fractures. The participants underwent daily gait assessments for several weeks using in-shoe IMU sensors. Eight gait parameters and each coefficient of variance were calculated. Some physical tests, such as hand grip strength and Timed Up and Go tests, were performed to check the baseline functional ability. RESULTS: The fracture group showed lower foot angles of dorsiflexion and plantarflexion in the swing phase. The receiver operating characteristic curve analyses revealed that a total foot movement angle (TFMA) < 99.0 degrees was the risk of subsequent fracture. CONCLUSIONS: We extracted the daily-life gait characteristics of patients with DRF using in-shoe IMU sensors. A lower foot angle in the swing phase, TFMA, may be associated with the risk of subsequent fractures, which may be effective in evaluating future fracture risk. Further studies to predict and prevent subsequent fractures from daily-life gait are warranted.

Lower-limb sagittal joint angles during gait can be predicted based on foot acceleration and angular velocity.

Background and purpose: Continuous monitoring of lower-limb movement may help in the early detection and control/reduction of diseases (such as the progression of orthopedic diseases) by applying suitable interventions. Therefore, it is invaluable to calculate the lower-limb movement (sagittal joint angles) while walking daily for continuous evaluation of such risks. Although cameras in a motion capture system are necessary for calculating lower-limb sagittal joint angles during gait, the method is unrealistic considering the setting is difficult to achieve in daily life. Therefore, the estimation of lower-limb sagittal joint angles during walking based on variables, which can be measured using wearable sensors (e.g., foot acceleration and angular velocity), is important. This study estimates the lower-limb sagittal joint angles during gait from the norms of foot acceleration and angular velocity using machine learning and validates the accuracy of the estimated joint angles with those obtained using a motion capture system. Methods: Healthy adults (n = 200) were asked to walk at a comfortable speed (10 trials), and their lower-limb sagittal joint angles, foot accelerations, and angular velocities were obtained. Using these variables, we established a feedforward neural network and estimated the lower-limb sagittal joint angles. Results: The average root mean squared errors of the lower-limb sagittal joint angles during gait ranged between 2.5°-7.0° (hip: 7.0°; knee: 4.0°; and ankle: 2.5°). Conclusion: These results show that we can estimate the lower-limb sagittal joint angles during gait using only the norms of foot acceleration and angular velocity, which can help calculate the lower-limb sagittal joint angles during daily walking.

Identification of characteristics of foot position and angle during swing phase in fallers using principal component analysis.

Identifying the characteristics of fallers is important for preventing falls because such events may reduce quality of life. It has been reported that several variables related to foot positions and angles during gait (e.g., sagittal foot angle and minimum toe clearance) differ between fallers and non-fallers. However, examining such representative discrete variables may not be sufficient to detect crucial information, which may be contained in the large portions of unanalyzed data. Therefore, we aimed to identify the comprehensive characteristics of foot position and angle during the swing phase of gait in non-fallers and fallers using principal component analysis (PCA). Thirty non-fallers and 30 fallers were recruited for this study. We performed PCA to reduce the dimensions of foot positions and angles during the swing phase and obtained principal component scores (PCSs) for each principal component vector (PCV), which were then compared between groups. The results revealed that the PCS of PCV3 in fallers was significantly larger than that in non-fallers (p = 0.003, Cohen's d = 0.80). We reconstructed waveforms of foot positions and angles during the swing phase using PCV3 and our main findings can be summarized as follows. Compared to non-fallers, fallers have a 1) low average foot position in the z-axis (i.e., height) during the initial swing phase 2) small average foot angle in the x-axis (i.e., rotation in the sagittal plane), during the initial swing phase, and 3) large variability in foot position in the y-axis (i.e., anterior/posterior position) during the initial swing phase. We can conclude that these are characteristics of gait related to fallers. Therefore, our findings may be beneficial for evaluating fall risk during gait using a device such as a shoe- or insole-embedded inertial measurement unit.

Differences in the coordination and its variability among foot joints during running in neutral foot and flatfoot.

Flatfoot is a well-known foot deformity, with a prevalence of 11.2%-29.0% among adults. Running injuries can occur in individuals with flatfoot; however, the underlying mechanism remains unknown. We investigated the coordination pattern and variability among foot joints while running by comparing participants with neutral foot and with flatfoot. Participants with neutral foot (n = 15) and flatfoot (n = 15) were asked to run at their preferred speed. Using the modified vector coding technique, the coupling angle between the foot joints, representing interjoint coordination, was calculated and categorized into four coordination patterns. The standard deviation of the coupling angle was computed to measure the coordination variability during the stance phase. There were no differences in the spatiotemporal parameters (speed, step length, and cadence) between the groups. In the sagittal rearfoot and sagittal midfoot coordination patterns, the flatfoot group showed a significantly greater proportion of anti-phase with proximal dominancy and a lower proportion of in-phase with proximal dominancy than the neutral foot group during early stance. Coordination variabilities between the sagittal rearfoot and sagittal midfoot (midstance), between the sagittal midfoot and sagittal forefoot (early stance), and between the frontal rearfoot and sagittal midfoot (midstance) were greater in the flatfoot group than in the neutral foot group. This may explain why those with flatfoot are likely to experience running injuries.

Estimation of lower-limb sagittal joint moments during gait using vertical ground reaction force.

Lower-limb sagittal joint moments during gait are important variables for evaluating the risk of disease progression, such as that of orthopedic diseases. Therefore, quantifying lower-limb sagittal joint moments during walking is important to continuously evaluate the risk of disease progression. A motion capture system and force plate are employed in the calculation of lower-limb sagittal joint moments during gait. However, they cannot be used during daily walking. Therefore, it is important to estimate these moments during walking from the vertical ground reaction force (vGRF), which can be measured using a wearable sensor, such as an insole device. Thus, this study aimed to estimate the lower-limb sagittal joint moments during gait using only the vGRF and confirmed its accuracy. This study included 188 healthy adults, and each participant walked at a comfortable speed (10 trials). We estimated the moments from the vGRF using a feedforward neural network. Our major findings are that our method can estimate lower-limb sagittal joint moments using the vGRF with accuracies of NRMSE¯ within 6.0-11.7% (NRMSEs¯ of the hip, knee, and ankle were 8.4, 11.7, and 6.0%, respectively). To the best of our knowledge, this study is the first to estimate lower-limb sagittal joint moments (including those of the hip, knee, and ankle joints) during gait using only the vGRF. Our method may be useful to estimate lower-limb sagittal joint moments during daily walking using only the vGRF, which can be measured by an insole device in the future.

Influence of sex and knee joint rotation on patellofemoral joint stress.

PURPOSE: Females are two times as likely to experience patellofemoral pain syndrome (PFPS) than males, however, the reason for this difference between sexes remains unclear. Patellofemoral joint (PFJ) stress is believed to contribute to PFPS alterations through knee joint rotation alignment, but the influence of knee joint rotation conditions on PFJ stress is unclear. We aimed to investigate the influence of sex and knee joint rotation alignment on PFJ stress. METHODS: Simulation ranges were set to knee joint flexion angles of 10-45° (common to both sexes) and extension moments of 0-240 Nm (males) and 0-220 Nm (females). The quadriceps force and effective lever arm length at the quadriceps muscle were determined as a function of the knee joint flexion angle and extension moment. The PFJ contact area, which is specific to sex, and knee joint rotation were calculated from cadaver data, and PFJ stress was estimated. RESULTS: In all knee joint rotation conditions, PFJ stress was higher in females than in males. Additionally, PFJ stress in males and females was the largest under neutral conditions compared with other rotation conditions. CONCLUSION: The results of the present study may be useful for understanding the underlying mechanisms contributing to the differences in PFPS in males and females.

Effects of step length and cadence on hip moment impulse in the frontal plane during the stance phase.

BACKGROUND: An excessive daily cumulative hip moment in the frontal plane (determined as the product of hip moment impulse in the frontal plane during the stance phase and mean number of steps per day) is a risk factor for the progression of hip osteoarthritis. Moreover, walking speed and step length decrease, whereas cadence increases in patients with hip osteoarthritis. However, the effects of step length and cadence on hip moment impulse in the frontal plane during the stance phase are not known. Therefore, this study aimed to examine the effects of step length and cadence on hip moment impulse in the frontal plane during the stance phase. METHODS: We used a public dataset (kinetic and kinematic data) of over-ground walking and selected 31 participants randomly from the full dataset of 57 participants. The selected participants walked at a self-selected speed and repeated the exercise 15 times. We analyzed the data for all 15 trials for each participant. Multiple regression analysis was performed with the hip moment impulse in the frontal plane during the stance phase as the dependent variable and step length and cadence as independent variables. RESULTS: The adjusted R 2 in this model was 0.71 (p < 0.001). The standardized partial regression coefficients of step length and cadence were 0.63 (t = 5.24; p < 0.001) and -0.60 (t =  - 4.58; p < 0.001), respectively. CONCLUSIONS: Our results suggest that low cadence, not short step length, increases the hip moment impulse in the frontal plane. Our findings help understand the gait pattern with low hip moment impulse in the frontal plane.

The Effects of Differences in the Morphologies of the Ulnar Collateral Ligament and Common Tendon of the Flexor-Pronator Muscles on Elbow Valgus Braking Function: A Simulation Study.

The anterior bundle (AB) and posterior bundle (PB) of the ulnar collateral ligament and the anterior common tendon (ACT) and posterior common tendon (PCT) of the flexor-pronator muscles have an independent form and an unclear form. The purpose of this study was to clarify the effect of differences in the morphologies of the AB, PB, ACT, and PCT on the elbow valgus braking function. This investigation examined three elbows. In the classification method, the AB, PB, ACT, and PCT with independent forms constituted Group I; the AB, ACT, and PCT with independent forms and the PB with an unclear form constituted Group II; the AB, PB, ACT, and PCT with unclear forms constituted Group III. The strains were calculated by simulation during elbow flexion at valgus at 0° and 10°. At 0° valgus, Group I and Group II showed similar AB and PCT strain patterns, but Group III was different. At 10° valgus, most ligaments and tendons were taut with increasing valgus angle. The average strain patterns of all ligaments and tendons were similar for the groups. The AB, PB, ACT, and PCT may cooperate with each other to contribute to valgus braking.

Shank and rearfoot coordination and its variability during running in flatfoot.

Flatfoot is a risk factor for patellofemoral pain syndrome (PFPS), and excessive rearfoot eversion occurring in flatfoot has been associated with the development and progression of PFPS; however, the mechanism remains unclear. This study aimed to investigate transverse shank and frontal rearfoot coordination patterns and variability when running with normal foot and flatfoot. Participants with normal foot (n = 13) and flatfoot (n = 13) were asked to run at their preferred speed. The coupling angle between the shank and rearfoot, representing intersegmental coordination, was calculated using the modified vector coding technique and categorized into four coordination patterns. Standard deviation of the coupling angle was computed as a measure of coordination variability during the stance phase. No differences in the characteristics and spatiotemporal parameters between groups were found, and all participants had rearfoot strike pattern. During midstance, the flatfoot group showed a significantly greater proportion of anti-phase with proximal (shank) dominancy than the normal foot group (p = 0.04, effect size = 0.88 [large]). Furthermore, flatfoot group showed a significantly greater in variability than the normal foot group (p = 0.03, effect size = 0.91 [large]). This study's results may help explain why flatfoot is likely to result in PFPS. However, the occurrence mechanism of running injuries like PFPS is multi- factorial. Since these results alone are not sufficient to explain the cause-effect relationship between flatfoot and injuries like PFPS, a prospective study including other factors such as patellofemoral joint stress would also be needed.

Morphological features of the lateral plantar ligament of the transverse metatarsal arch.

The aim of this study was to elucidate the morphological characteristics of the lateral Lisfranc ligament in a large sample. This investigation examined 100 legs from 50 cadavers. Each of the lower limbs was dissected to identify the plantar aspect of the transverse metatarsal arch, and morphological characteristics of the lateral plantar ligament were assessed, including the length, width, and thickness of the fiber bundles. The majority of plantar ligaments originated from the base of M5 and the plantar aspect of the lateral cuneiform (LC). The lateral plantar ligament could be classified into three types: Type I, a band-like fiber bundle originating from the base of M5 to the LC (41%); Type II, originating from the base of M5 and the plantar aspect of LC and mostly connected the blending the fiber bundles of the tibialis posterior (TP) and long plantar ligament (LPL) (21%); and Type III, with no ligaments originating from the base of M5 and plantar aspect of the LC (38%). The morphological characteristics of Type I lateral plantar ligament were as follows: length, 31.8 ± 3.7 mm; width, 2.3 ± 1.0 mm; and thickness, 0.2 ± 0.3 mm. The morphology of the lateral plantar ligament showed variation, originating from the base of M5 and the plantar aspect of LC most commonly, but this was not the case in 38% of limbs. The findings suggest that the lateral plantar ligament might play a role in the transverse tarsal arch, indicating a cooperative mechanism with the TP and LPL.

Differences in rearfoot, midfoot, and forefoot kinematics of normal foot and flatfoot during running.

Flatfoot is a common foot deformity, which could contribute to running injuries such as medial tibial stress syndrome. Intrafoot kinematics of flatfoot during walking have often been documented using multisegment foot models. However, the intrafoot kinematics of flatfoot during running remains unclear, despite the possible relationship between flatfoot and running injuries. We aimed to clarify rearfoot, midfoot, and forefoot kinematics when running in participants with normal foot and flatfoot. Participants with the normal foot (n = 14) and flatfoot (n = 14) were asked to runover-ground at their preferred speed. Three-dimensional kinematics of the rearfoot, midfoot, and forefoot during running were calculated based on the Rizzoli foot model. A two-sample t-test of statistical parametric mapping was performed to determine differences between normal foot and flatfoot in time histories of intrafoot kinematics during running. No differences were found between groups in characteristics and spatiotemporal parameters. In the frontal rearfoot angle, a significantly increased eversion from 24% to 100% (p < .001) was observed in the flatfoot compared to the normal foot. At the midfoot angle, a significantly increased eversion from 0% to 4% (p < .049) and 21% to 100% (p < .001) was observed in the flatfoot compared to the normal foot. At the forefoot angle, a significantly increased inversion from 6% to 17% (p < .047) was observed in the flatfoot compared to the normal foot. These findings may be useful to explain why flatfoot could contribute to running injuries such as medial tibial stress syndrome.

A new method for estimating three-dimensional movement of the patella using a surface mapping method and computed tomography.

INTRODUCTION: A previous study reported a method called the 2D-3D registration technique to examine three-dimensional movement of the patella. However, that method requires a biplane fluoroscopy system. In the present study, the aim was to establish a new method (CT-based surface mapping method) to estimate three-dimensional positions and angles of the patella with a motion capture system and CT. METHODS: In Study 1, the most appropriate parameters for the CT-based surface mapping method (i.e., target edge length, threshold of thickness of the soft tissue, and minimum distance between markers) were explored and determined. In Study 2, three-dimensional movement (i.e., positions and angles) of the patella using the CT-based surface mapping method and the most appropriate parameters were determined, and they were compared with the true positions and angles obtained by CT. RESULTS: The results of Study 1 showed that the most appropriate conditions were as follows: (1) target edge length, 3 mm; (2) threshold of thickness of the soft tissue, 0-20 mm; and (3) minimum distance between markers, 10 mm. The results of Study 2 showed that the errors of the positions and angles were less than approximately 10 mm and 10° at most, respectively (both supine and sitting positions). CONCLUSION: The CT-based surface mapping method may be useful for a future study to clarify differences in three-dimensional movements of the patella between patients with patellar tendinitis and healthy subjects.

Effect of Gender and Load Conditions on Foot Arch Height Index and Flexibility in Japanese Youths.

Arch height index (AHI) and arch height flexibility (AHF) are useful methods for evaluating foot structure. Although foot structure may be linked to intrinsic factors such as gender and load conditions, information on AHI and AHF in consideration of these factors is lacking. This study aimed to examine the effect of gender and load conditions on AHI and AHF. One hundred Japanese youths (50 males, 50 females) were recruited in this study. Arch height and truncated foot length were measured with an AHI measurement system. AHI was calculated for each load condition using truncated foot length and arch height. AHF was defined as the change in arch height from 10% to 50% of weightbearing load, and from 10% to 90% of weightbearing load. To satisfy the assumption of independence, only measurements from the right foot were analyzed. A gender × load condition interaction was found in AHI. AHI in all load conditions showed significant differences between the genders (p < .001), and AHI of female participants was significantly less than that of male participants (p < .001). In contrast, no significant gender × load condition interaction was noted in AHF, and only the main effect of the load condition was found (p < .001). In AHI use, the effects of gender and load conditions must be considered, and AHF may be used considering only changes in load conditions. The results of this study provide useful information regarding which normative values of AHI and AHF should be used.

Anatomical variations in the insertion of the peroneus longus tendon.

BACKGROUND: This study aimed to inspect anatomical variations in the insertion of the peroneus longus tendon (PLT) using a large sample of cadavers. METHODS: In total, 104 legs from 52 Japanese cadavers were used. The PLT was identified behind the lateral malleolus and carefully followed up to its insertion in the foot. All insertion slips of the PLT were located and documented. RESULTS: Mainly, the PLT was inserted to the base of the first metatarsal (1MT) in all 104 ft. Attachment to the medial cuneiform was present in 20.2%, and the first dorsal interossei was present in 36.5%. The anterior frenular ligament was observed in 31.7%, and attachment to the flexor digiti minimi brevis and opponens digiti minimi was present in 31.7%. The posterior frenular ligament was observed in 5.8%. An additional band was observed in 3.9%, and the adductor hallucis consisting of a caput obliquum was present in 3.9%. No statistically significant differences in the PLT were observed between genders or laterality (right vs. left). CONCLUSIONS: These findings suggest that the main function of the PLT is resisting the varus force on the 1MT; however, as the PLT has various attachment sites, it may also be involved in the stabilizing action of the longitudinal and transverse arches. Therefore, these variations and functions appear to be associated with a difficult diagnosis at the first clinical evaluation.

Algorithm to compute muscle excitation patterns that accurately track kinematics using a hybrid of numerical integration and optimization.

Forward dynamic simulation is used to examine the causal relationships between muscle excitation patterns and human movement. The computed muscle control (CMC) algorithm computes a set of muscle excitations for a movement using proportional-derivative control. However, errors between experimental and simulated kinematics may cause rapid movements. Herein, we propose a novel algorithm, i.e., hybrid computed muscle control (HCMC), which uses a hybrid of numerical integration and optimization to compute muscle excitation patterns that accurately track kinematics, even for rapid movements. We compared the muscle excitation patterns and accuracies of the kinematics simulated by HCMC and CMC using synthetic and experimental data. Two simple musculoskeletal models were used. The synthetic data were generated for three repetitive movements from the rest position to the flexed position (the hip, knee, and ankle underwent 10°, 20°, and 10° plantar flexion, respectively) and back to the rest position for various times. Experimental data were obtained for a subject running at 220 steps/min. The maximum errors in all kinematics calculated using the HCMC algorithm were extremely lower than those calculated using CMC algorithm (HCMC: 0.04-0.07° [synthetic data] and 0.00-0.03° [experimental data]; CMC: 1.04-2.41° [synthetic data] and 0.48-2.50 [experimental data]). For rapid movements, muscle excitations estimated using HCMC occurred early and without delay than those estimated using CMC. The HCMC algorithm can provide muscle excitation patterns that accurately track kinematics and may be useful for perturbation studies using forward dynamic simulation of joints characterized by a low range of motion during rapid movements.

A mathematical modelling study investigating the influence of knee joint flexion angle and extension moment on patellofemoral joint reaction force and stress.

BACKGROUND: Patellofemoral pain (PFP) is the most common orthopaedic condition among runners. Individuals with PFP exhibit greater patellofemoral joint (PFJ) reaction force and stress when compared with pain-free controls. However, it is not clear whether PFJ reaction force and stress are the highest (or lowest) when knee joint flexion angle and extension moment are in which combinations. We aimed to investigate the influence of knee joint flexion angle and extension moment on PFJ reaction force and stress. METHODS: A PFJ sagittal model was used to quantify PFJ reaction force and stress. Based on the public dataset of the previous study, peak knee joint flexion angle and extension moment at various running speeds was calculated. Based on the calculated peak value, simulation ranges were set to knee joint flexion angle of 10-45° and extension moment of 0-240 Nm. The quadriceps force, effective lever arm length at quadriceps muscle, and PFJ contact area were determined as a function of the knee joint flexion angle and extension moment, and finally PFJ forces and stress were estimated. RESULTS: PFJ reaction force increased as the knee flexion angle and extension moment increased. Although PFJ stress also increased as the knee extension moment increased, it was at the highest and lowest at 10° and about 30° knee joint flexion angles, respectively. CONCLUSIONS: Incorporating knee flexion posture (approximately 30°) during running may help in reducing PFJ stress, which would be useful in the prevention of pain and act as an optimal treatment program for PFP.